Analysis of Strain Fields in Silicon Nanocrystals 
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Strain has a crucial effect on the optical and electronic properties of nanostructures. We calculate 
the atomistic strain distribution in silicon nanocrystals up to a diameter of 3.2 nm embedded in 
an amorphous silicon dioxide matrix. A seemingly conflicting picture arises when the strain field 
is expressed in terms of bond lengths versus volumetric strain. The strain profile in either case 
shows uniform behavior in the core, however it becomes nonuniform within 2-3 A distance to the 
nanocrystal surface: tensile for bond lengths whereas compressive for volumetric strain. We reconsile 
their coexistence by an atomistic strain analysis. 



The low dimensional forms of silicon embedded in 
silica have strong potential as an optical material, [lj 
Such heterogeneous structures inherently introduce the 
strain as a degree of freedom for optimizing their opto- 
electronic properties. It was realized earlier that strain 
can be utilized to improve the carrier mobility in bulk 
silicon based structures. [2j This trend has been rapidly 
transcribed to lower dimesional structures, starting with 
two-dimensional silicon structures. @ Recently for silicon 
nanowires, there have been a number of attempts to tailor 
their optical properties through manipulating strain. [3,0 
Futhcrmore, recent studies have revealed that the strain 
can become the major factor restricting the crystalliza- 
tion of the nanolayers. [f| 0] It depends on several factors, 
most important of which are the lattice mismatch be- 
tween the constituents, size of the NCs, and the growth 
conditions, such as the details of the growth procedure. 
In summary, for improving the optical and electronic 
properties of nanocrystals (NCs), the strain engineering 
has become an effective tool to be exploited. [H, yuj El A 
critical challenge in this regard is to analyze the strain 
state of the Si NCs embedded in silica. 

The close relations between strain and optical or elec- 
tronic properties in Si NCs have very recently become 
the center of attention. [13, [H, G3 There still remains 
much to be done in order to understand strain in nanos- 
tructures at the atomistic level. As pioneered by Tsu et 
al. Raman spectroscopy can be an effective experimental 
tool for determining the strain state of the Si NCs.fl4j 
Specifically, recent Raman studies reported that the Si 
NCs may be under a thermal residual strain and this can 
be reduced by overall annealing at high temperatures [f| 
or by local laser annealing. [l5j] Due to small density dif- 
ference between Si NC and the surrounding a-SiC>2, a 
limited information can be gathered about its structure 
using transmission electron microscopy (TEM) or even 
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high resolution TEM techniques. [16] Especially, molecu- 
lar dynamics simulations with realistic interaction poten- 
tials present an opportunity, by providing more detailed 
critical information then the best imaging techniques cur- 
rently available and clarify the analysis of experimental 
results. Along this direction, previously [l7| we focused 
on Si-Si bond length distribution and reported that Si-Si 
bond lengths are stretched upto 3% just below the sur- 
face of Si NCs embedded in amorphous Si02 which has 
also been very recently confirmed. [1 S| 

In this Letter, we analyze the volumetric and bond 
length strain distributions in Si NCs, in particular 
demonstrate that both compressive volumetric strain and 
tensile bond length strain coexist within the same Si 
NC. We accomplish this by performing trajectory anal- 
ysis on model samples (with ca. 5000 atoms) simulated 
via molecular dynamics u sing a reliable and accurate as 
well as reactive force field. |19l| T he simulation details are 
similar to our previous work,[l7[ except the way we con- 
struct the Si NC in glass matrix. Instead of deleting 
all glass atoms within a predetermined radius, we re- 
move the glass atoms after rigorously defining the sur- 
face of the nanocrystal through the Delaunay triangula- 
tion method. [20j In this way, we have constructed NCs 
embedded in glass matrix with diameters ranging from 
2.2 nm to 3.2 nm without introducing built-in strain to 
the system. In this diameter range we observe similar 
trends in strain, volumetric strain, and bond length dis- 
tribution etc., therefore, we present only the figures of 
the system for a typical NC of radius 2.6 nm. 

In the language of geometry, strain is defined through 
an affine transformation that maps the undeformed state 
to deformed state, which is called deformation gradient. 
Several methods to derive discrete form of deformation 
gradient from atomic positions are reported. [U [Hj] 
In the method proposed by Pryor et al., the atomistic 
strain tensor is derived from local transformation ma- 
trix that transforms nearest neighbors of a certain atom 
from its undeformed state to the deformed one. From the 
MD simulations, using positions of NC atoms, we first 
extract each atoms displacement vector from its unde- 



Distance to the Surface (A) 

FIG. 1: (Color online) Variation of, Si-Si bond lengths 
(squares), hydrostatic strain (diamonds), and the volumet- 
ric strain (triangles) as a function of distance to nanocrystal 
surface (see text). Dashed, dotted and solid lines are guides 
to the eye for the respective data set. The data for 2.6 nm di- 
ameter NC is used. Inset: Other NC diameters ranging from 
2.2 nm to 3.2 nm are also shown. 
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FIG. 2: (Color online) Dependence of solid angle subtended 
by tetrahedron face to the angle between tetrahedron face 
and nanocrystal surface. Illustration of solid angle subtended 
by tetrahedron face (top left inset) and the angle between 
tetrahedron face and NC surface (bottom right inset). 



formed site which is determined by positioning an ideal 
tetrahedron to the local environment. Using these dis- 
placement vectors, we construct deformation matrix and 
derive the atomistic strain tensor from this local defor- 
mation tensor. (23j The first invariant of strain tensor cor- 
responds to the hydrostatic strain. [HI As an alternative 
measure to hydrostatic strain, we calculate volumetric 
strain by considering volume change of a tetrahedron 
from its undeformed counterpart. A third measure as we 
have used in our previous report, [l7[ is the bond length 
strains. 

To verify our results we have calculated strain distri- 
bution in NC region for all mentioned measures. We have 
plotted all three of them in Fig. [TJ The results of volu- 
metric strain are very close to hydrostatic strain which 
is the trace of strain tensor calculated with aferomen- 
tioned technique. [23j] In these results, we observe a net 
compressive behavior of strain just under the surface and 
a uniform tensile strain of about 1% at the core of NC. 
Si-Si bonds are stretched by about 1% in the core region 
in agreement with the hydrostatic and volumetric strains, 
however, just under the surface, Si-Si bonds are stretched 
up to 3% where hydrostatic and volumetric strain results 
indicate compressive strain state. The bond-stretch in Si- 
Si bonds due to oxidation has been shown earlier by us 
using molecular dynamics simulations 17] which was also 
confirmed by other approaches, [l^ Occurrence of com- 
pressive volumetric strain and stretched bond lengths in 
the same outer region may initially seem contradicting. 
However, stretching of bonds does not imply that the 
system is under tensile hydrostatic strain as well. Con- 
sider a tetrahedron formed by a Si atom and its four 
Si neighbors (A, B, C, D) as shown in upper inset of 
Fig. [21 In the ideal case, the solid angle (O) subtended by 



each triangular face of this tetrahedron should be equal 
to 180°. Under a uniform deformation, bond lengths 
will also be stretched, while the solid angles remain un- 
changed. However, under a nonuniform deformation, the 
change in three solid angles causes a decrease in the vol- 
ume of the tetrahedron while increasing or preserving the 
bond lengths. Hence, a combination of stretched bond 
lengths with deformed solid angles may end up with an 
overall reduction of the volume of the tetrahedron. This 
explains the coexistance of compressive volumetric strain 
and stretched bond lengths at the region just below the 
surface of NCs. 

To better visualize the nature of the deformation of 
the Si NCs, we consider the orientational variation of the 
solid angles of the tetrahcdral planes. As illustrated in 
the lower inset of Fig. [21 the two important directions are 
the unit normal (ns) of the tetrahedron face subtending 
the solid angle under consideration, and the local out- 
ward surface normal (tt-nc) of the NC. It is clearly seen 
from Fig. that solid angles subtended by tetrahedra 
faces oriented outward to the NC surface are increased 
up to 220°, whereas those facing inward to the NC core 
are decreased down to 160°. This dependence is a clear 
evidence of how oxidation affects strain distribution close 
to the interface. 

To further quantify the atomistic strain in the highly 
critical region within 3 A distance to the interface, we 
classify the average bond length and hydrostatic strain 
behaviors into three categories. Figure [3] displays the 
percentage as well as the bonding details of each cate- 
gory. In top- left, we illustrate most common type with a 
share of 53.0% which is responsible for the opposite be- 
havior in Fig. [T] where average bond lengths of center Si 
atoms to its four nearest neighbors are stretched but net 
atomistic strain at this atom is compressive. In this case 



3 




Hydrostatic Strain = -2.7% A V Hydrostatic Strain = 2.0% 
Av. Bond Length = 5.2% * Av. Bond Length = 2.6% 




Hydrostatic Strain = -3.6% 
Av. Bond Length = -1.2% 

FIG. 3: (Color online) Illustrations of oxidation effects on 
strain in three categories with their percentage of occurrences: 
Si-Si bonds are stretched and system is under compressive 
strain (upper left). Si-Si bonds are stretched and system is 
under tensile strain (upper right). Si-Si bonds are shortened 
and the Si atom at the center is under compressive strain (bot- 
tom). Large spheres (gold) and small spheres (red) represents 
Si and O atoms, respectively. 

solid angles facing toward the oxide region is increased 
to 270° due to oxygen bonds of Si neighbors. Although 
these oxygen bonds stretched Si-Si bonds to 2.41 A, net 
strain on center Si atom is -2.7%. In the top-right part 
of the Fig. [3] we illustrate second most often case with 
a percentage of 42.0%, where average bond lengths and 
atomistic strain show similar behavior; bond lengths are 



stretched and net hydrostatic strain is tensile. In this 
case oxidation somewhat uniformly deforms the bonds 
so that solid angles are still around 180° which is the 
value for the unstrained case. Finally, as shown at bot- 
tom of Fig. [3l a very small percentage of atoms (5.0%) 
in the region beneath the surface have shortened bond 
lengths and compressive atomistic strain. 

In summary, we study the strain state of Si NCs in sil- 
ica matrix with diameters in 2 to 3.2 nm. The structure is 
assumed to be free from any thermal built-in strains. The 
core region of the NC is observed to be under a uniform 
1% tensile strain, where both bond length and volumet- 
ric strain measures are in agreement. However, towards 
the NC interface, while the Si-Si bonds become more 
stretched, the hydrostatic strain changes in the compres- 
sive direction. In the interpretation of the indirect strain 
measurements eg. from spectroscopy, this dual character 
needs to be taken into consideration. We explain these 
two behaviors using the solid angle deformation of the 
tetrahedral-bonded Si atoms, and demonstrate that it is 
ultimately caused by the oxygen atoms at the interface. 
An equally important finding is that the overall strain 
profile within the Si NCs is quite nonuniform. As very 
recently emphasized, within the context of centrosym- 
metric materials, like silicon, such strain gradients locally 
break the inversion symmetry and may lead to profound 
physical consequences. [25] 
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